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ABSTRACT 

We present 21 cm H I observations of four Hickson Compact Groups with evidence for a substantial 
intragroup medium using the Robert C. Byrd Green Bank Telescope (GBTJ^- By mapping H I emis¬ 
sion in a region of 25'x25' (140-650 kpc) surrounding each HCG, these observations provide better 
estimates of H I masses. In particular, we detected 65% more H I than that detected in the Karl 
G. Jansky Very Large Array (VLA) imaging of HGG 92. We also identify if the diffuse gas has the 
same spatial distribution as the high-surface brightness (HSB) H I features detected in the VLA maps 
of these groups by comparing the H I strengths between the observed and modeled masses based on 
VLA maps. We found that the H I observed with the GBT to have a similar spatial distribution as 
the HSB structures in HCGs 31 and 68. Conversely, the observed H I distributions in HCGs 44 and 
92 were extended and showed significant offsets from the modeled masses. Most of the faint gas in 
HCG 44 lies to the Northeast-Southwest region and in HCG 92 lies in the Northwest region of their 
respective groups. The spatial and dynamical similarities between the total (faint-l-HSB) and the 
HSB HI indicate that the faint gas is of tidal origin. We found that the gas will survive ionization by 
the cosmic UV background and the escaping ionizing photons from the star forming regions and stay 
primarily neutral for at least 500 Myrs. 

Subject headings: 


1. INTRODUCTION 

Hickson Compact Groups (HGGs) are dense concen¬ 
tration of 4-10 galaxies that have densities si milar to 
galax y clusters, but are relatively isolated (iHicksonl 
Il982fl . Their low veloci ty dis persions (AU ~200 km 
s“^ ; see iHickson et al.l (ll992tH and very high galaxy 
densities favor extreme tidal interactions. Such in¬ 
teractions alter the morphology and distribution of 
atomic gas in their host galaxies and have been 


studied by extensively ('' 

A^illiams van Gorkoml 

1995; 

Verdes-Montenegro et al. 

20011: IWilliams et al.l 

2002; 

Verdes-Montenegro et al. 

20021 120051 iDurbala et al. 

20081 Ide Mello et al.1 l2012al). The atomic gas in some 


of these groups exists in complex tidal structures and 
even spatial offsets from their host gal axies. The groups 
are also found to be deficient in H U (jWilliams fc Rood! 


Electronic address: sanch@pha.jhu.edu 

^The National Radio Astronomy Observatory is a facility of the Na¬ 
tional Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 

^ This is referred to as “H I -deficiency” and is de¬ 
fined as Defni = log[M(HI)p^e(i] - log[M(HI)obs] (see 


ll987HHnchtmeieHll99111. 

IVerdes-Montenegro et al.l (l200lL VMOl hereafter) sug¬ 
gested that the nature and distribution of H I in the 
tidal structures can be used as a proxy to estimate the 
level of interaction and hence the ‘phase of evolution’ of 
the group. They proposed that the groups in their early 
stages of evolution would have galaxies with individual 
H I envelopes. As the group evolves, the tidal interac¬ 
tions would increase and result in severe tidal stripping. 
As a result, most of the H I in such systems should ex¬ 
ist as tidal debris. In the final stages of evolution the 
stripped H I gas from multiple galaxies gets stirred-in 
together. This blends the individual clumps and facil¬ 
itates the formation of large H I envelopes around the 
groups. It is also likely that during the process of in¬ 
teractions part of the neutral gas goes through a phase 
change via ionization, thus making the evolved groups 
deficient in H I. 

In a study bv iBorthakur et ^ (|20I0l BIO hereafter), 
the authors compared H I data obtained with the Karl 

IHavnes &: Giovanellil[T984l . for predicted HI masses). 
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Fig. 1.— Digital Sky Survey r-ba nd images of H CGs 31, 44, 68, and 92 in greyscale with 21 cm VLA Hi maps overlaid as black contours. 
The member galaxies as defined by IHicksM ll 198211 are labeled for extended groups (HCG 44 and 68). The GBT beam size (FWHM=9.1') 
for our pointings are overplotted in grey. 


G. Jansky Very Large Array (VLA) to that from the 
Robert C. Byrd Green Bank Telescope (GBT). They 
found that the groups have a significant fraction of H I in 
a faint and diffuse phase that was missed by the VLA. 
They also found that the fraction of H I in the diffuse 
component increases with the evolutionary phase of the 
group. They advocated that with increasing interaction 
time as marked by the evolutionary phase of the group, 
the tidally-stripped H I structures disperses more and 
more into a faint H I medium, which finally may even 
get ionized. This was identified as a contributing factor 
for the observed H I deficiency. 

The presence of substantial amount of neutral gas in 
the intragroup medium can have long-lasting influence on 
the evolution of the group. The high neutral column den¬ 


sities and pressures resulting from tidal dynamics have 
been suggested as a p otential cause for ind ucing star for¬ 
mation in tidal debris (iMullan et al.ll201^ . Several stud¬ 
ies have found large numbers of young stellar clusters be¬ 
ing formed in the intragroup me d ium (iHu nsberge r et al.l 
119961 : iTorres-Flores et ^ 120091 : iGallagher et al.l I2010L 
and citations therein). For example the intergalactic 
star-forming regions in Stephan’s Quintet appear to be 
devoid of older (> 10® yrs) stellar populations and 
their SED resembles that of dusty star-forming regions 
in the galactic disk s than dwarf star-forming galaxies 
l|Boauien et al.l[20T^ . Although, the diffuse gas detected 
by BIO is unlikely to be the sites of star formation, nev¬ 
ertheless the diffuse gas may eventually feed star forming 
regions in the group. On the other hand, the starburst- 
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Fig. 2.— GET 21cm HI spectra for each of the four HCGs in our sample. The spectra have been plotted in a grid format similar to the 
5x5 pointing scheme on the sky. Pointings are separated by 4' from their neighbors resulting in a total coverage of 25'x25'. The center 
of the map coincides with the center of the group. The data presented in this paper have been smoothed to 10 km s“^ . The position of 
the pointings are marked at the top left corner and their observed HI masses at the top right corner. We have listed the 5 (t HI masses for 
spectra where no HI was detected. The modeled HI masses are shown in color. They are color-coded to show the agreement/disagreement 
between the observed and modeled masses in green/red. Pointings where the limiting H I masses are consistent with the predicted masses 
are shown in blue. The systemic velocity of the member galaxies are marked, (a) 25 spectra towards HCG 31. The depression centered at 
4200 km s~^ is an artifact left over after clearing RFI at that frequency. The group size as defined bv IHicksonI 11198211 and re-defined by 
IVerdes-Montenegro et ^ II2005I 1 places the entire group well within the central pointing. The redshifts of the member galaxies place all of 
them within the velocity coverage of the HI profile. At the redshift of the group, each pointing is separated by 66 kpc. The broad feature 
between 3170-3600 km s“^ has no associated optical galaxy. It may be a RFI feature although not clear from the current observations. 
Further observations are required to confirm it. 


driven winds can significantly impact the morphology of 
the optical and the H I structures and may even pro¬ 
duce extreme diffe rences between the H I and optical 
tidal morphologies (|Hibbard et al.l[2000ll . The winds will 
be most effective against the low column density diffuse 
gas and can ionize it thus enhancing the H I deficiency. 
Therefore, it is crucial to understand the distribution and 
conditions in the diffuse gas to predict its fate and its role 
in the overall evolution of the group. 

Motivated by these studies, we present our pilot study 
of the extent of the H I in four HGGs using the GET. 
Based on results by BIO, we infer that all the four groups 
have a significant fraction of their H I in the faint diffuse 
extended phase thus making them ideal candidates for 
studying the extended nature of H I. We start by present¬ 
ing our sample followed by details of the observations, 
and data analysis in Section 2. The results and a discus¬ 
sion of their implications are presented in Section 3. Fi¬ 
nally, we conclude our findings in section 4 and comment 
on the scope of a future study. The cosmological pa¬ 
rameters used in this study are Hq = 70 km s“^ Mpc“^, 
Om = 0.3, and Oa = 0.7. 


2. OBSERVATION 
2.1. Sample 

Our sample for the pilot study consists of four groups: 
HGGs 31, 44, 68, and 92^. These four groups were part of 
the original sample of BIO and were found to have a sub¬ 
stantial fraction of their H I in the faint diffuse extended 
phase. All of the groups have at least 4 spectroscopically 
confirmed members and are within 100 Mpc of the Milky 
Way. An advantage of this sample is that these groups 
have ext ensive H I data includin g single dish observa - 
tions by IWilliams fc RoodI (|1987ll . iHuchtmeie^ (Il997ll . 
and BIO as we ll as 21 c m H I i maging with the VLA by 
IWillianis et al.l (TlOOll. l2002h : IVerdes-Montenegro et al.l 
( 200,^ : lTang et M.l (|2008| 1 and Westerbork S ynthesis Ra¬ 
dio Telescope (WSRT) bv ISerra et al.l (|2013f) . The prop¬ 
erties of the sample including their redshifts, distances, 
group sizes, our GBT map centers, physical scales corre¬ 
sponding to 1', and references to published VLA H I im¬ 
ages are presented in Table 1. The optical images show- 

^ HGG 92 is known as Stephan’s Quintet. Our observations do 
not include the foreground galaxy HGG 92A (aslo known as NGC 
7320) at a redshift of 0.002622 (v = 786 km s”! ) 
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Fig. 2.— (b) Spectra from 23 pointings around HCG 44. Pointings are separated by 4' or 23 kpc at the redshift of the group. We 
were unable to obtain spectra for two of the pointings due to time constraints during observations. Features associated with the individual 
galaxies are identified and are marked A, C, D, or J corresponding to 44A, 44C, 44D, and SDSS J1012. 


ing VLA H I contours are presented in Figure [TJ 

The H I distribution of these groups can be classified 
into two distinct morphologies. The H I in HCGs 31 and 
92 exists as complex tidal structures indicative of intense 
tidal interactions in these groups. VLA imaging showed 
that a large fraction of the H I in these two groups is not 
associated with individual galaxies (sometimes even spa¬ 
tially offset from the galaxies), but exists as tidal debris 
in the group environment. Therefore, hereafter we refer 
to them as tidal debris-dominated groups. The H I in 
HCG 44 and 68 is mostly concentrated as interstellar 
medium (ISM) of individual member galaxies and hence 
they will be referred to as ISM-dominated groups for the 
remainder of the paper. It is worth noting that even the 
ISM-dominated groups are dehcient in H I , although the 
H I in the VLA images does not show any major tidal 
structures. This is probably indicative of a quick phase 
change of the neutral material in the intragroup medium 
(IGrM). The evolutionary model classifies HCG 31 to be 
in phase 2 of evolution and the remaining three groups, 
HCGs 44, 68, and 92, in phase 3. 

2.2. Observations and Data Reduction 

We obtained 21cm H I maps with the GET under pro¬ 
gram GBT07A-093. The motivation behind these obser¬ 
vations was to probe the H I distribution in these selected 
HCGs beyond the region surveyed by BIO and to explore 
the spatial distribution of the diffuse gas detected by BIO. 
Our objective was to ascertain if the H I in the diffuse 
phase has similar spatial distribution as the high-surface 


brightness (HSB) gas detected in the VLA maps. We 
conducted the observations in pointing map mode with 
25 individual pointings on a 5x5 grid with a spacing of 4' 
between pointings. The total surveyed area was 25'x25'. 
The positions of the centers of the maps are presented 
in Table [1] and the positions of the GBT pointings are 
shown in Figure [1] as gray circles. We employed stan¬ 
dard GBT observing routines for the mapping and spent 
a fifth of the time on the reference position. The refer¬ 
ence position was chosen on a blank part of the sky at 
an offset of 1 degree from the edge of the map in both 
Right Ascension (RA) and Declination. 

To ensure positional accuracy we performed local 
pointing corrections using the observing procedure Au- 
toPeak. Hence, the error in pointing should be less than 
GBT’s global pointing accuracy of 3". This is less that 
0.6% of the full width at half maximum (FWHM) of 9.1' 
for the GBT beam at 1.4 GHz. The flux was calibrated 
using L-band calibrators 3G 48 (16.5 Jy), and 3C 286 
(15.0 Jy) and the antenna gain was derived for each ses¬ 
sion individually. We state a calibration error of ~ 10% 
(Scr) in all the measurements presented in the paper. 

The spectroscopic data were acquired in dual polariza¬ 
tion L-band and cover the frequency range between 1.15 
and 1.73 Gigahertz (GHz). Using 9-level sampling and 
two IF settings, each with 12.5 MHz total bandwidth, 
the spectrometer gives 8192 channels at 1.5 kHz (0.3 km 
s“^) resolution, covering a total velocity range of 2500 km 
s“^. Data were recorded every 10 seconds in order to iso¬ 
late radio frequency interference (RFI) and to remove the 
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Fig. 2. — (c) Spectra showing the Hi distribution around HCG 68. The pointings are separated by 4' or 38 kpc at the redshift of the 
group. The only galaxy with HI in this group is HCG 68C. It shows an unperturbed doubled-horned profile and has a large extended disk 
of in the VLA image. The peak of the profile centers at position B4, which is closest pointing to the position of 68C. 


corrupted integrations during data processing. Integra¬ 
tion times per pointing, the range of velocities over which 
H I was detected, and range of errors in H I masses are 
provided in Table [T] 

Data for two more groups (HGGs 15 and 100) were 
also obtained under the same program. Unfortunately, 
they were severely corrupted by extensive RFI. There¬ 
fore, these two groups will not be discussed any further. 

The data were analyzed using the interactive soft¬ 
ware package GBTIDL, specifically designed for reduc¬ 
tion and analysis of spectral line data obtained with the 
GBT. Gomposite spectra were obtained by combining 
each RFI-free integration using GBTIDL routines and 
the final spectra were fitted with a low-order polynomial 
(second or third order in most cases) for estimating the 
baseline. The good bandwidth stability has allowed us 
to fit a lower order polynomial and thus made the data 
sensitive to broad faint spectral features. In addition we 
chose a large range in velocity to fit for the baseline and 
hence we expect to detect individual features that are 
narrower than ^1500 km s“^ . 

The reduced spectra are presented in Figure [2] in grid 
format. The position of the grid is marked on the top 
left corner. The labels A-E mark the decrease in RA 
and the numbers 1-5 mark the decrease in Declination. 
Adjacent pointings are not independent due to the large 
beam size of the GBT (FWHM=9.1'). Table [T] presents 
distances between pointings with respect to the central 
pointing i.e. the center of the map. We provide figures 
showing the GBT spectra overlaid on the VLA maps as 


additional images in the journal version. 

3. RESULTS AND DISGUSSION 

Atomic gas was detected in all four groups, although 
their properties vary widely. Here, we interpret out data 
in conjunction with the GBT spectra for these groups 
from our earlier study (BIO) and VLA imaging of these 
groups (see Tabled] for references). In our earlier study, 
we found that the GBT spectra had similar spectral 
shapes as those of the VLA for most of the groups, 
although with higher amplitudes. The only additional 
features detected in the GBT spectra that were not 
seen in the VLA data were H I “wings” that connected 
distinct velocity components with fainter emission (e.g. 
HGG 92). 

In this study, we focus on the variation of spectral pro¬ 
files as a function of distance from the center of the group. 
Our GBT map covers four times the group size (see Ta¬ 
ble dj) for HGGs 31 and 92 (where extensive H I tidal 
structures were observed in VLA images) and double 
the group size for HGG 68. HCG 44 is the only group 
where the mapped area is not significantly larger than 
the group. For this group, we refer the reader to an ex¬ 
tended mterferomfitricmap of the region surrounding the 
group (|Serra et al.ll2013ll . 

3.1. Spectral Signature of Atomic Hydrogen 

The two ISM-dominated groups, HCGs 44 and 68, 
show distinct H I features at the systemic velocities of the 
member galaxies. For example, HCG 68 shows a double¬ 
horned profile at the systemic redshift of HCG 68C (NGC 
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Fig. 2.— (d) Spectra showing the Hi distribution around HCG 92. The pointings are separated by 4' or 107 kpc at the redshift of the 
group. The redshift of the member galaxies are marked The peak flux was observed at position C2 and the observed Hi mass at C2 was 
used to normalize the modeled masses. 


5350) peaking at the pointing closest to its spatial lo¬ 
cation. Similarly, the spectra for HCG 44 show dis¬ 
tinct features at the optical redshifts of individual mem¬ 
bers HCG 44A (NGC 3190), 44C (NGC _ 3185 ), 44D 
(NGC 3187) and SDSS J1017 (iSerra et aHl2?fill . The 
spectral regions corresponding to the member galaxies 
are labeled in Figure [2] 

The tidal debris-dominated groups, HCGs 31 and 92, 
do not exhibit distinct H I features corresponding to 
the spatial locations and redshifts of the member galax¬ 
ies. HCG 31 exhibits a single component H I feature 
(AVflux>o A 350 km s“^ and ct = 57 km s“^) encom¬ 
passing the redshift range of all the group members. 
H I emission in HGG 31 peaks at the center of the 
map. The shape of the profile remains the same with 
a slight shift of less than 50 km s“^ in the peak as we 
move from north to south. There is no change in veloc¬ 
ity of the peak from east to west. Towards the edge of 
the map in pointings Al and A2, another strong emis¬ 
sion feature was detected between 3170-3600 km s“^ . 
There is no corresponding galaxy found in the vicinity 
of this detection in the Sloan Digital Sky Survey (SDSS) 
image. This feature is dynamically not connected to the 
H I associated with HGG 31 and is separated by more 
than 300 km s“^. In addition, it was only detected in 
pointings Al and A2, and not in Bl, B2, and A3. This 
suggests that the spatial location of this H I structure is 
to the Northeast of HGG 31, outside our 25'x25' map. It 
is possible that this feature is not an astronomical source, 
but was generated by RFI. However, we do not see a de¬ 


tectable variation in its strength between integrations as 
well as the two polarizations. Observations from a sepa¬ 
rate epoch are required to confirm this detection. 

The spectra towards HGG 92 show three features. The 
strongest feature covers a velocity range of 6350-6750 km 
s“^ and peaks at position G2. This matches the velocity 
range of some of its member galaxies: HGG 92E and 92D, 
and extends as far out to the redshift of 92G. On the low 
velocity end, the H I feature has a large H I “wing” first 
detected in our previous study (BIO) that can be seen 
extending out to the second H I feature in the higher 
signal to noise spectrum. The second feature lies between 
5950-6150 km s“^ . It peaks at the redshift of galaxy 
HGG 92F and in positions G2, G3, B2, and B3. The third 
feature is seen between velocities 5600-5800 km s“^ , thus 
covering the redshift of HGG 92B and peaks in position 
B3. 

All the groups with the exception of HGG 68 show 
blending between various H I components. Most of 
the H I in HGG 44 exists in stable disk-like structures 
(double-horned profiles). The proximity of the systemic 
velocity between the member galaxies within the groups 
does result in blending in some of the pointings. How¬ 
ever, we found the spatial variation among spectral fea¬ 
tures to be strongest for the ISM-dominated groups. The 
H I kinematics in the tidal-debris dominated groups are 
much more complex. There are no double-horn disk¬ 
like signatures in HGGs 31 and 92 (except for 92F). In 
fact, most of the H I in both the groups exists as a sin¬ 
gle peaked profile covering the redshift of all/some the 
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Fig. 3.— Integrated Hi spectra of all the pointings around the four groups are shown in black (see § 3.2). For comparison we show the 
spectrum from the pointing with most HI flux in red. For HCG 44, the spectrum from the central pointing is shown in blue and is of 
similar strength as the pointing of maximum flux. The labels on the top right hand corner specify the group, its total HI mass, and the 
HI mass associated with the pointing containing the most flux. The peak flux is 98% of the total flux observed in HCGs 31 and 68. This 
implies that most of the gas is in the central region and gas in the outskirts is just 2% of the total HI in these groups. In HCG 92 the gas 
is more spread out and the peak is 81% of the total HI whereas in HCG 44 the peak flux is just 54% of the total flux. The broad feature 
in the spectra of HCG 31 between 3170-3600 km s“^ is not associated with any optical galaxy in the SDSS images and are certainly not 
asscociated with the group. Hence that feature was left out of the mass estimated printed on the top right-hand corner. 


member galaxies. The tidally dominated groups shows 
stronger kinematic blending between various components 
in all pointings. Both these observations are in agree¬ 
ment with the evolutionary model which predicts that 
in the final stage of evolution (phaseSb), H I , if sur¬ 
viving outside galaxies, would form a single cloud en- 
veloping the entire group with a single velocity profile 
(jVerdes-Montenesfro et al]l200lh . 

3.2. Mass of the Atomic Hydrogen 

The atomic gas mass for each of the pointings is pre¬ 
sented in Table [T] The masses are also printed on the 
top right corner of each of the panels in Figured These 
measurements were made using the relationship M(HI) = 
2.36x10^D^(SAV) Mq, where D is the luminosity dis¬ 
tance to the galaxy in megaparsecs (Mpc) (column 4 in 
Table[T|) and SAV is the velocity integrated H I flux den¬ 
sity in Jy km s“^ . We also estimated the total FI I flux 
spread over our 5x5 grid by adding spectra from each 
of the pointings and then normalizing it by a factor of 
(tt In2)/Z^ = 5.86, where Dfwhm = 9.1' is the 

FWHM of the GBT beam and I = 4.0' is the length of 
each grid on our map. The resulting spectra are shown 
in Figure [3] and the associated FI I masses are presented 
in the last row of Table [TJ 


The four groups cover a wide range of H I masses. 
Within the surveyed area, the tidal debris-dominated 
groups FICG 31 and 92 have 1.86 x 10^° and 2.45 x 
10^° Mq of H I mass respectively whereas the ISM- 
dominated groups, FIGG 44 and 68, contain significantly 
less atomic gas - 1.82 x 10® and 8.32 x 10® Mq respec¬ 
tively. These integrated FI I masses are close to the to¬ 
tal H I associated with the groups except in the case of 
FIGG 44 where the FI I extends beyond the GBT mapped 
area. The total H I mass of the FIGG 44 as reported by 
ISerra et all (|2013L WSRT imaging) was 2.2 x 10® Mq. 

On the other hand, the H I mass detected in the 
VLA imaging toward s HCG92 is merely 1.47 xlO^® Mq 
(IWilliams et al.l[200M including Anon2 and corrected for 
cosmological parameters). This indicates that our GBT 
measurements found 65% more FI I in this group than 
that observed in the VLA image. Similarly, for FIGG 
68 the GBT detected FI I mass is 10% large r than that 
detec ted in the VLA map, 8.32x10® Mq (|Tang et al.l 
l2008l) . Both GBT and VLA confirm that the FI I in 
FIGG 68 is contained in the double-horn profile associ¬ 
ated with galaxy 68G. Also, in FIGG 31 the GBT detected 
60% more FI I than that dete cted in CnB configuration 
VLA image of 1.126x10^® Mq (|Verdes-Montenegro et all 
120051 cosmological parameters corrected). The new GBT 
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Fig. 4.— Models illustrating the variation in observed flux (and 
mass) as a function of distance between the source and the pointing 
center when observed with the GBT at 1.4 GHz. The GBT beam 
is a Gaussian of FWHM of 9.1' and hence is sensitive to flux at 
>3% level for distances upto 12'. Three source geometries are 
considered - point source (red), a square source of side 5' (green), 
and a circular source of diameter 9' (blue). The sensitivity of 
the GBT at large offsers from the beam center is higher for more 
extended targets as indicated by the blue curve representing the 
circular source with the larger area among the three geometries 
considered here. We note that these simulations are for symmeteric 
sources. Large asymmetries in the source geometry can result in 
significant deviations in the flux sensitivity profile. 

measurements also match that of single di sh measure- 
ments by BIO and iWilliams fc RoodI (|l987f) with a dif¬ 
ference of less than 5%. However, the DnC configuration 
image of HCG 31 was reported to have a larger total 
H I mass (IVerdes- Montenegro et all 120051) . Part of the 
reason for the difference could be the result of the seri¬ 
ous solar interference that effected the DnC configuration 
data extensively. 

For comparison, we plot the pointing with the maxi¬ 
mum flux (see Figure [2]) in red in Figure |31 In all the 
groups the total flux is larger than the maximum seen 
in an individual pointing. The difference in flux be¬ 
tween the total and the peak measurement ((M(HI)totai — 
M(HI)max)/M(HI)totai) Is about 2% for HCGs 31, and 68; 
19% for HCG 92; and 46% for HCG 44. In all the groups 
the gas on the edges (M(HI)edges) is dynamically similar 
to the central pointing. This is indicative of large spa¬ 
tial coverage of the GBT beam (well beyond its FWHM 
of 9.1'). The total H I spectra of HCGs 92 and 44 also 
show additional H I features that were not detected in 
the pointing with maximum flux. These features add 
substantial amount of H I to the total atomic gas mass 
of the groups. Interestingly, the third emission feature 
in HCG 92 observed between 5700-5800 km s“^ peaks 
sharply in the total flux spectrum and not in any of the 
individual spectra. This is suggestive of a large spatial 
extent of the gas associated with this feature. 

It is noteworthy that while tidal debris-dominated 
groups contain multiple galaxies, their spectra, when 
observed with a large beam, is indistinguishable from 
individual systems. This is a valid concern for the 
blind H I surveys that do not distinguish between gas 
associated with individual galaxies and gas in galaxy- 
groups. For example, when observed with a beam of 
FWHM « 15', the compactness of such groups would 


make it difficult to identify them as complex H I tidal 
structures. However, such H I -rich groups are very 
rare both in terms of galaxy distribution as well as 
in terms of their H I masses. The two tidal debris- 
dominated groups would exist at th e higher-mass end o f 
the H I mass function obtained by IZwaan et al.l (120051 1 
for the nearby universe bas ed on the H I Parkes All Sky 
Survey (iMever et al.l l2004 HIPASS;)^. In the case of 
HGG 92, only 2% of the extragalactic H I detections 
have mass greater than that of this group. On the other 
hand, the ISM-dominated groups lie close to the peak of 
the H I mass function, i.e. 10®'® Mq. The above com¬ 
parison illustrates the gas richness of these systems. 

3.3. Spatial Extent of the Atomic Hydrogen in HCGs 

Our data consist of 25 spatially distinct pointings. 
However, they are not independent. This is because 
the pointings are separated by 4' whereas the GBT 
beam is a two-dimensi onal Gaussian of FWHM = 9.1' 
(|Boothrovd et al.l 1201 111 . As a result, the GBT beam is 
sensitive to sources at ^ 5% levels at even 11.5' away 
from the pointing center. Figure |4] illustrates the effect 
of the GBT beam shape on its sensitivity as a function 
of distance between the source and the pointing cen¬ 
ter. To further illustrate the effects of source geometry 
on the sensitivity of the GBT, we also considered three 
source geometries - point source, source in the shape of 
a square of side 5', and a circular source of diameter 
9.1'. The sampling in these simulations matches that 
from our observations. These simulations imply that 
none of our pointings are completely independent of the 
central pointing. This also suggests that the strength of 
the flux observed at any position is dependent on the 
shape/morphology of the source. 

A comparison of our data with the models shown in 
Figure m are presented in Figure O The plots show the 
radial distribution of H I as a function of distance from 
the central pointing in the units of physical distance. The 
observed masses are shown in red circles. The modeled 
H I masses are for two test cases: a point source and 
a circular source of diameter 9.1'are shown as a green 
dashed line and a brown dotted dashed line, respectively. 
Unresolved systems are close to the green dashed line 
and barely resolved systems are sightly above the brown 
dotted-dashed line. The H I structures in HGG44 and 
HGG 68 are resolved whereas in HGG 31 they are un¬ 
resolved. Also, H I structures in HGG 92 are partly 
resolved. We also show the models based on the H I dis¬ 
tribution seen in the VLA maps (discussed later) as blue 
diamonds. 

To disentangle the effects of the beam, we imaged the 
data using the GBT Mapping Pipeline. However, the 
sampling of data, i.e. 25 pointings, was found to be sig¬ 
nificantly lower for studying distribution and morphology 
of the gas. We also tried to extracted spatial informa¬ 
tion from the data using the Richardson-Lucy algorithm 
to deconvolve the GBT beam pattern. Unfortunately, 
the results were found to be of low spatial resolution 
and consistent with those found with the GBT Mapping 
Pipeline. 

In the absence of high resolution maps of the groups 
with the GBT, what can we learn about the H I distri- 

^ average gridded HIPASS beam is Ril5.5b 


















Faint Atomic Gas in HCGs 


9 


10.5 
10.0 
X 9.5 


M 9.0 


8.5 

8.0 




GBT mass HCG31 
Simulated masses: 
VLA HI map 
Point source 
Circular source ■ 


"i 


50 100 150 200 

Radius (kpc) 


250 


9.0<i. 


S 8.5 


O 

J 8.0 


7.5 


■ii 


GBT mass HCG44 
Simulated masses: - 
VLA HI map ; 
Point source 
■ ■ Circular source - 

o 

H : 


20 40 60 80 

Radius (kpc) 


100 


GBT mass HCG68 i 
Simulated masses: i 
VLA HI map : 
Point source 
Circular source - 



50 100 

Radius (kpc) 


O 


10.5 

10.0 

9.5 
9.0 

8.5 


k 






<>s 


GBT mass HCG92 
Simulated masses: 
VLA HI map 7 
Point source 
Circular source - 


V <:> 


100 200 300 400 

Radius (kpc) 


Fig. 5.— Radial distributions of Hi in the four groups The solid red circles represent the observed H I mass and the blue open diamonds 
represent the modeled masses. The ranges in the observed and modeled masses are shown as vertical bars. Each of these values was obtained 
by averaging 4-8 pointings that are equidistant from the center of the map. The exception to this is the value at zero which represents the 
central pointing. The models for a point source and a circular source as shown in Figure |4]are overplotted in green and brown respectively. 
It is worth noting that not all the pointings that were averaged are independent, although the modeled masses takes that into account. 
The model data were normalized at the central pointing (C3) corresponding to zero in the absissa. The data from HCG 31 matches that 
of a point source whereas the other sources are slightly resolved. The range of values in the observed and the modeled data match well for 
HCG 31 and 68 whereas HCG 44 and 92 show large differences in the spread. 


bution in these groups? It turns out that although we 
cannot image the H I , we can still ask whether the GBT 
data are consistent with the spatial distribution of HSB 
H I seen by the VLA. For instance, if the excess gas that 
was missed by the VLA but detected by the GBT (BIO), 
has a similar spatial distribution then we expect to see 
the GBT data agree with the VLA data. Alternatively, 
if the excess gas has a very different spatial distribution, 
then the GBT and the VLA maps should be inconsistent. 

We simulated the expected masses for each of our GBT 
pointings using the VLA H I distribution as a template 
for the source structure. The hypothesis that we are test¬ 
ing here is that the total H I (faint -f-HSB) distribution 
traced by the GBT has a similar spatial distribution to 
the HSB features detected by the VLA. Since the GBT 
is sensitive to more diffuse gas in addition to the HSB 
features detected by the VLA, we normalized the mod¬ 
eled masses to one of the pointings. In this case, we used 
the pointing with the maximum observed mass. This 
is an arbitrary choice and normalization can be done 
at any pointing. The normalization process assumes a 


“constant” value of the ratio total-to-HSB H I for the 
entire group. So in this case the “constant” for normal¬ 
ization was chosen to be the ratio at the pointing of peak 
H I flux/mass. 

A good match between the observed and modeled 
masses would mean that the ratio is constant across 
the mapped region. In other words, this would imply 
that the total (faint-l-HSB) H I distribution traced by 
the GBT has a similar spatial distribution as the HSB 
features detected by the VLA. Gonversely, any deviation 
would indicate that the assumed ratio is not constant 
spatially, thus implying a different spatial distribution of 
the faint gas. It is worth noting that since the modeled 
masses were normalized, the deviation can be both pos¬ 
itive and negative. A deficit in observed mass does not 
correspond to the missing H I in the GBT observations. 
Instead it simply corresponds to a smaller total-to-HSB 
H I ratio (also equivalent to smaller faint-to-HSB H I ra¬ 
tio) than that of the pointing with maximum flux. 

Figure [2] shows the observed masses in black at the top 
right-hand corner and the modeled masses in color. De- 
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viations between the observed and the modeled masses 
that exceed more than 0.1 dex (= ±20%, which is twice 
the uncertainty in the measurements) were considered 
significant and to be beyond the noise in our measure¬ 
ments. The modeled masses for the pointings, where 
they are in agreement with the observed masses, are 
printed in green (i.e. |Mmod — Mobsl < 0.1 dex). Point¬ 
ings where the modeled masses differ from the observed 
masses by more than 0.1 dex are indicated in red (i.e. 
iMmod — Mobsl > 0.1 dex). Pointings where the mod¬ 
eled masses are smaller than the limiting H 1 masses (i.e. 
Mmod ± Miim) are printed in blue. 

3.3.1. HCG 31 

The observed and the modeled H 1 distributions are 
consistent for most of the pointings for HCG 31. This 
indicates that most of the excess gas seen in the GBT 
spectra is associated with HSB features and has a similar 
spatial distribution. Hence, we conclude that most of the 
H I is centrally located and is consistent with being asso¬ 
ciated with the ISM/tidal structures of member galaxies. 
The dispersion in the observed masses as a function of 
radial distance and that of the modeled masses suggest 
that the VLA and the GBT spatial distribution are sim¬ 
ilar. In fact, the observed distribution is consistent with 
a point source (relative to the GBT beam) as seen in 
Figure [5l We have also analyzed the observed H I dis¬ 
tribution by deconvolving the GBT beam from the data 
using the Richardson-Lucy algorithm and found consis¬ 
tent results. 

Pointing El shows the maximum difference between 
the observed and modeled masses corresponding to a 
mass excess of M(Hl)obs — M(Hl)sim = 7 x 10® Mq. 
This points to the presence of an H I structure at the 
edge of the surveyed region in the Southwest direction 
such that none of the other pointings have detected it 
at >20% levels. We also observed small discrepancies 
between the modeled and observed masses in pointings 
Bl, D5, and E4. Unlike El, the observed masses in these 
pointings were smaller than the modeled masses. This 
indicates that the faint H 1 distribution is more skewed 
in the Northeast to Southwest direction than the HSB 
features. 

3.3.2. HCG 44 

The H 1 in this group shows distinct spatial vari¬ 
ations corresponding to the positions of the member 
galaxies. All the H 1 features peak at the position of 
their host galaxies and gradually fall off with distance. 
The strongest H I feature is associated with HGG 44D 
(NGC 3187) and peaks at position B4. This pointing was 
used to normalize our modeled masses. Since the VLA 
map did not cover channels beyond 1800 km s“^ , the 
spectral feature between 1900-2000 km s“^ associated 
with galaxy SDSS J1012 (marked as ‘J’ on the spectra) 
was not covered in the VLA maps. The total H I mass, 
including the feature corresponding to SDSS J1012, is 
printed in black at the top right corner in Figure 3. For 
pointings B5, G4, G5, D4, and D5, we also present the 
H 1 masses for gas at n < 1800 km s“^ in grey in paren¬ 
thesis. This is the most appropriate quantity to compare 
with the modeled fluxes estimated from the VLA map. 

Most of the pointings in the Northwest to South¬ 
east direction show significant excesses from the mod¬ 


eled masses. This ratio of total to HSB H I increases 
rapidly with distance from B4. The ratio, (M(HI)gbt — 
M(Hl)sii„)/M(Hl)sini, peaks at pointings E5 (^ 700%) 
and Cl (~ 600%). All the pointings towards the North¬ 
east are in agreement. This suggests that most of the 
faint gas lies towards the South of B4. This is not sur¬ 
prising as two of the Hl-rich neighbors of 44D are South 
and Southwest of its position. 

Figure [5] shows the radial distribution of H I in 
HCG 44. The normalization was done at the central 
pointing for ease in comparison. The profile clearly de¬ 
viates from a point source or even circular source filling 
the central 9.1'. H I associated with this group extended. 
Interestingly, the H I radial profile of this group shows 
significant differences in the dispersion between the ob¬ 
served and modeled masses. 

3.3.3. HCG 68 

The H I detected in this group comes from the double¬ 
horn profiled centered at the redshift of galaxy HCG 68C 
(NGC 5350). The spatial position of H I is consistent 
with originating in the H I disk of this blue spiral galaxy 
HCG 68G. The H I peaks at position B4 that is closest 
to the position of galaxy HGG 68G. The good match 
between observed and modeled masses thus suggests that 
the total H I distribution is almost identical to the HSB 
features observed by the VLA. 

The H I profile and distribution are consistent with 
a large unperturbed H I disk as imaged b y the VLA 
and s hown in Figure [1] (also see Figure 25 of iTang et al.l 
l2008ll . The profile is highly symmetric and shows no sign 
of dynamical tidal deformation. The H I diameter of the 
disk is ^ 9' from the VLA H I maps. At the redshift of 
the galaxy, this implies a disk size of ^86 kpc in diameter. 
It is rare for such an extended disk to show no sign of tidal 
interaction in a group environment. Therefore, galaxy 
HGG 68G is most likely a recent addition to the group 
and could still be infalling. 

3.3.4. HCG 92 

The relative strengths of the three features associated 
with the H I in HGG 92 show variations with positions 
suggesting the presence of distinct features at different 
positions and velocities. The strongest feature peaks at 
position C2 suggesting it is offset from the center of the 
group by more than 2'. The velocity range of this fea¬ 
tures matches that with Arc-N, Arc-S, and NW-HV that 
were identified in th e VLA imaging (see Figure 8 from 
I Williams et al.l [200^ . The positions of these structures 
are consistent with the offsets seen in our GBT data. 

Figure [2] shows the distribution of modeled masses and 
its comparison to observed masses. The modeled fluxes 
are normalized at pointing G2. We found that the distri¬ 
bution is consistent in the center of the group, however it 
varies significantly at the outskirts namely at pointings 
Gl, Bl, A2, A3, B4, G4, D4, D3, D2, Dl, and E2. While 
the first four pointings show an excess in the observed 
H I masses, the remaining seven pointings show a deficit. 
This represents a displacement of the total H I distri¬ 
bution from the HSB features towards the North of the 
group. This is also consistent with the fact that the 
high-velocity H I wing between 6300-6530 km s“^ , first 
observed by BIO, is primarily localized between point¬ 
ings B2, B3, G2, and G3. This suggests that the gas 
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associated with the H I “wing” is located close to the 
intersection between these pointings. 

3.4. Nature and Origin of the Diffuse Extended Gas in 

HCGs 

The spatial distribution of the diffuse gas is similar to 
HSB features in HGGs 68 and 31. This implies that 
the diffuse gas is in faint extensions of the HSB gas. 
On the other hand, HGGs 92 and 44 shows that most 
of the diffuse gas is concentrated in a localized region 
and has an overall distribution different from that of the 
HSB features. The discrepancy (excess & deficits) be¬ 
tween observed and the modeled masses suggests that 
the diffuse gas has a different spatial distribution than 
the HSB structures. The spectral profile of the total 
H I being similar to HSB features suggests that the dif¬ 
fuse gas may be associated with older tidal structures. 
For example, HGG 44 shows dynamically disturbed HSB 
H I structures near the pointings with the maximum dis¬ 
crepancy at Cl and E5. Similarly, the location of the gas 
associated with the high-velocity “wings” (flux between 
6300-6530 km s“^ seen only in the GET data) is con¬ 
sistent with the Arc-N, which exhibits not only complex 
morphological and dynamical structure, but is also spa¬ 
tially displaced from any of the member galaxies, thus 
indicating a tidally disrupted structure. 

The existence of faint diffuse structures associated 
with tidal interactions in multi-galaxy interacting sys¬ 
tems is well known. One such example is the extended 
tidal structures s een in the nearby galax y-group, M81 
(|Yun et al.l 11994 IChvnoweth et al'l l2008ll . Strong dy¬ 
namical disruption in the spectral profiles of the cold 
gas in HCGs 44 and 92 indicate their origin to be related 
to the transformation of the ISM of individual galaxies 
into the intragroup medium. 

3.5. Fate of the Gas and its Role in the Evolution of the 

Group 

In this sub-section, we discuss the survival timescales 
of diffuse structures against photoionization from the ul¬ 
traviolet (UV) background and ionizing radiation pro¬ 
duced by young star forming regions in the groups. The 
H I in the ISM remains immune to ionization by the back¬ 
ground UV radiation owing to the high densities inside 
their host galaxies. However, after being stripped, the 
ISM clouds are unable to maintain their internal pressure 
and hence expand due to their internal velocity disper¬ 
sion. As a result the density and the column density of 
the clouds drop. This can be expressed as a function of 
time, t, and the internal velocity dispersion, v, as 

n = ni LI/{Li + vt)^ and N = Ni L^/{Li + vt)^ (1) 

where, Lj, n^, and are the initial sizes, densities, 
and column densities. Consequently, the tidally stripped 
H I structures expand with time and become more and 
more susceptible to ionization from the cosmic UV back¬ 
ground. 

They can also be ionized by extreme ultraviolet radi¬ 
ation that may be escaping from the star forming re¬ 
gions (if any) in their vicinity. Tidal interactions are 
known to trigger the formation of young stars. At times 
some of the starforming regions may be hosted at the 
tips of extended tidal tails and resemble dwarf galax¬ 
ies. Such structures have been appropriately termed 



Log [Density (cm’^)] 


Fig. 6.— Variation of column density for an Hi cloud as a function 
of density that is being irradiated the cosmic UV background. The 
model was generated using photoionization code CLOUDY for a 
cloud with a hydrogen column density of 2x10^^ cm~^, shown as 
the green dots. The blue and the red dots refer to neutral and 
ionized hydrogen as a function of the density of the cloud. As 
the density drops, the cloud becomes increasingly susceptible to 
ionization by the cosmic UV background. The ionization state 
of the cloud changes from primarily neutral to primarily ionized 
between a density of 10~^ cm~^ and 3 x 10“^ cm~^. 


as tidal-dwarf galaxies (TDG: IZ wickvl [ l956t iSchweizerl 
Il978t iMirabel et al.l 119921 : iDuc &: Mirabell Il999ll . Like 
other interacting galaxies, HCGs have been found to host 
multiple star-forming TDGs a nd individual star clus¬ 


ters outside the parent galaxies (iHunsberger et al.l 


Ig]esias;Haramo^yilhie3j200]£jGaJlaghfir^^_aljJ^U^^ 

Xu et al.ll2005HTorres-Flores et alJl2009tlGallagher et al l 


1996 


2001 


2010l:lde Mello et al.ll2012bHKonstantoDOulos et al.ll2013t 

Eigenthaler et al.ll2015HFedotov et al.ll2015L and others'). 
The young stars produce large amounts of ionizing pho¬ 
tons a nd given the right conditions (such as str ong stellar 
winds I Vogt et ^120131 : iBorthakur et al.ir2014ll a fraction 
of the ionizing flux may escape the molecular cloud co¬ 
coon where the young stars are embedded. 


3.5.1. Ionization by the Cosmic UV Background 

In order to investigate the lifetime of such tidally 
stripped neutral gas structures, let us consider photoion¬ 
ization of the neutral gas by the UV background. Ob¬ 
servations suggest that atomic gas at column densities 
greater than 2 x 10^® cm~^ is immune to ionization from 
the U V background (|Corbelli et al.l 119891 : Ivan GorkomI 
Il991h . This result is based on deep 21 cm H I mapping 
of the disks of normal galaxies that reveal a sharp cutoff 
in the disks at this H I column density. Motivated by 
this argument, we can consider the time it would take a 
stripped gas cloud to reach this column density. This will 
be its survival time. For example, a stripped ISM struc¬ 
ture of initial size I kpc with density of 1 cm“^ expanding 
at 20km s“^ would take «500 Myrs to reach this column 
density. At that time, the density will drop by more than 
three orders of magnitude to 7.5 x 10“^ cm“^. If it con¬ 
tinues to expand beyond this time then the structure will 
be susceptible to ionization. Hence we can conclude that 
such a H I structure will survive for at least 500 Myrs. 
On the other hand, the density of the h ot IGrM is ex- 
pected to be 4 X 10“^ — 9 x 10“^ cm“^ ([Freeland et al.l 
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Fig. 7. — Variation of column density of neutral and ionized hydrogen as a function of ionization parameter, Usb defined in Eq. 2), for 
a cloud. The Usb encompasses cloud properties such as its density and distance from the radiation source. We present two photionzation 
models that were generated using code CLOUDY. The left panel shows the ionization for a cloud of total hydrogen column density of 
2x10^^ cm~^ and the right panel for a column density of 2x10^^ cm~^ respectively. The total column density is kept constant (shown 
in green) while the variation of the neutral and ionized components are shown in blue and red respectively. As the ionization parameter 
increases so is the ability of the radiation field to ionize the cloud. As expected the partial ionization is achieve at lower ionization parameter 
for low column densities as can be seen from left panel as opposed to the right panel. The ionizing source for both the runs is a young 
stellar population (continuous starburst) of 4 Myrs of age. 


I2008I [Rasmussen et al.l[200l) . Therefore, pressure equi¬ 
librium requires that tidally stripped cooler ISM clouds 
will cease to expand much before reaching densities com¬ 
parable to that of the IGrM. Hence, the clouds would 
continue to survive for a much longer period of time. 

To explore the fate of the diffuse gas further, we 
performed detaile d photoioniza t ion m odeling using the 
code CLOUDY (|Ferland et al.l Il998[ l. where we con¬ 
firmed that clouds with neutral column density greater 
than 2x10^® cm“^ survive photoionization from the UV 
background. These simulations used the latest correc¬ 
tions to the UV background i.e., ionization parameter, 
U=<1)hm/(c uh), where $ hm = 2750 photons cin s“^ 
(similar to that used by iBorthakur et al.ll2013ll. b ased 
on the publishe d value s bv iHaardt fc Madad l]2012l l and 
ICiroiix fc Shulll (jl997ll . Even for clouds at N(H) = 
2 X 10^® cm“^, the ionization fraction due to photoioniza¬ 
tion would depend on the density of the cloud (see Fig¬ 
ure El). The ionization fraction changes drastically for 
clouds with densities between 10“^'® and 10“^ ° cm“^. 
Clouds with densities greater than 7.5 x 10“^ cm“^ have 
H I/Htotai fraction of more than 50%. In conclusion, 
as long as the column densities do not drop below this 
value, the diffuse H I structures will be able to survive 
ionization by the UV background. 

It is worth noting that the timescale calculations are 
dependent on the assumption for initial column density 
of the gas. Here we have chosen the initial values as the 
average density and column density expected in the spi¬ 
ral arm of a normal galaxy. However, the higher column 
density knots will be much more immune to ionization 
and may produce TDC or star-forming clumps. Such 
clumps are also likely to be detected in the VLA imaging 
due to their higher column densities. These clouds may 
survive indefinitely if they are not disrupted by other 


forces such as winds from young star-forming regions or 
further tidal disruption. Nevertheless, the large extent 
of the diffuse gas and their fairly long survival times sug¬ 
gest that the TDCs and the diffuse gas may have the 
same origin. For example, a back of the envelope calcula¬ 
tion shows that if a tidal arm is pulled out of the galaxy 
at the rotational speed of ^ 200 km s“^ for 500 Myrs 
then it will be displaced by 1 00 kp c from its host galaxy. 
Incidentally Ide Mello et al.l (120081 1 found TDCs in the 
H I envelope of HCG 100 at about 100 kpc away from 
the member galaxies. This indicate that the origin of 
TDCs and diffuse gas in some of the groups may be via 
the same process(es). Furthermore, the diffuse faint gas 
may eventually get accreted by the any of the galaxies, 
TDC, or larger star-forming clumps and help in sustain¬ 
ing star-formation in the group. 

3.5.2. Ionization by Young Stellar Populations 

The two tidal-debris dominated groups (HCGs 31 & 
92) show strong evidence for the pre sence of young 
stellar populations of 4-10 Myrs old (l}Ci et al.l 1200,4 
iTorres-Flores et akl 120091 : I Gallagher et al.ll2010ll . Here, 
we explore the possibility of photoionization of the low 
column density H I by the escaping ionizing radiation 
from such young star-forming regions in the vicinity of 
the H I in the intragroup medium. The effect of ioniz¬ 
ing radiation on the clouds can be characterized with the 
ionization parameter, Usb, which relates the properties 
of the starburst (SFR and escape fraction, /) and the 
H I cloud ( distance from the star burst, d, and density, 
nn) as Isee iTumlinson et al.ll201lh 


/ SFR N 

/ / \ /20 kpcN 

VIO Mq yr“i/ 

VO.Ol/ \ d J 


The higher the ionization parameter, the greater is the 
potency of the radiation field to ionize the gas. One of 
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Fig. 8. — Variation of the peak ionization parameter that a 
cloud can withstand while remaining primarily neutral, UsBmaxias 
a function of the total column density the cloud. UsBmax is pro¬ 
portional to the column density of the gas over more than 2 orders 
of magnitude. As the column density of the clouds increase, so 
does their ability to withstand a stronger radiation field. 


the important parameters influencing the ionization pa¬ 
rameter is the distance of the cloud from the radiation 
source. The closer the clouds are to the starburst regions, 
the higher their susceptibility for ionization. Therefore, 
given the properties of the cloud, we can estimate the 
minimum distance, Dhi from the starburst where H I will 
be able to withstand ionization and remain primarily 
neutral i.e. HI/{HI + HII) > 50%. Goming back to 
the case of an expanding tidal structure with initial col¬ 
umn density of 3 x 10^^ cm“^ and size 1 kpc, Dhi can 
be expressed in terms of the column density as 


log( 


=-1.198-0.5 log UsBu 


-0.75 


log( 


Nh 


3 X 10^^ cm“2 


+ 0.5 


SFR 


/ 


10 M0yr“i 0.01^ 


(3) 

where, UsBmax represents the largest value of the ion¬ 
ization parameter that the cloud can withstand while 
remaining primarily neutral. 

In order to explore the relationship between 
UsB, UsBmax, and the column density of the cloud, we 
modeled starburst driven photoionization of clouds us¬ 
ing the photoionization code, GLOUDY. The modeling 
assumes the geometry of the clouds to be in the form 
of a rectangular slab illuminated by a radiation field. 
A spectrum of a y oung stellar populati on of 4 Myrs 
(from Starburst99; iLeitherer et al.l 1201011 was used as 
the radiation field. At this age the feedback is strongest 
when the Wolf-Rayet and the O stars dominate the 
ionizing radiation and the stellar feedback^. Models 
were run for clouds with varying total hydrogen column 
density from 10^® to 10^^'® cm“^. Figure [7] shows 
the variation of column density of neutral and ionized 
hydrogen for two of the models. The UsBmax is identified 
by interpolating the column density between consecutive 
runs where neutral H I column crosses the 50% mark. 
Figure [8] shows the variation of UsBmax as a function of 
the column density of the clouds. As expected, as the 
column density of the clouds increase, so is their ability 
to withstand stronger radiation fields. 


However, similar results were obtained for a slightly older pop¬ 
ulation of 10 Myrs. 


Here we apply the results from our photoionization 
analysis to determine the minimum distance of the clouds 
to remain neutral for two of our groups, HCGs 31 and 
92, that are forming stars at a rate of 10 Mf:^ yr~^ 
(jGallagher et al.l[20in i and 6.7 Mq yr“^ (|Xu et al.ll2005ll 
respectively. A ssuming an escape fraction of ionizing 
photons as 1% (|Grimes et al.ll200^ . the minimum dis¬ 
tance or Dhi for clouds of column density 3 x 10^^ cm“^ 
is ^360 and ^^290 pc from the star forming regions of 
HCG 31 and 92 respectively. However, the same for 
clouds of column density 2 x 10^® cm“^ are ^340 and 
^280 kpc for HCG 31 and 92 respectively. 

HCG 31 is a very compact group and our data rule out 
the scenario of the faint gas being more than 300 kpc 
from the starforming regions. Then how can we under¬ 
stand the presence of faint low column density gas in 
this system? The likely answer is that faint gas is as¬ 
sociated with the HSB (high column density) gas that 
shields it from the starburst generated radiation field. 
This scenario is supported by the fact that the diffuse 
gas found by BIO has the same spectral distribution as 
the HSB features. On the other hand, in HCG 92, Dhi 
for a clouds of column density 5 x 10^® cm“^ is 77 kpc, 
which is smaller than the optical size of the group. In 
fact, our data on HCG 92 show that the diffuse gas in the 
H I wing (6300-6500 km s“^ ) is most likely associated 
with Arc-N, which is 80-100 kpc from the starforming 
regions. 

We note that besides photoionization, the starburst- 
driven winds might also lead to heating of the gas via 
shock ionization. However, the shock ionized gas (if 
any) will cool back r apidly as cooling time is small 
(jBorthakur et al.ll2013f l and it is unlikely that the shock 
ionization would contribute to the long term heat¬ 
ing/destruction of the neutral clouds. Detailed hydro- 
dynamical simulations will further help in probing this 
phenomenon with accuracy. In addition, analytical argu¬ 
ments by BIO indicate that the gas will also survive the 
conductive heating from the hot IGrM. The contribution 
of ram pressure stripping and heating f rom the hot X-ray 
intra group medium likely to be small (jR,asmussen et ahl 
[200l . In fact, even in most X-ray bright groups (ex - 
cept for the most massive ones'). IPesiardins et al.1 (120131 1 
found the X-ray emission is concentrated mostly around 
individual galaxies and does not arise from the group 
medium. 

Therefore, it is possible for the diffuse gas to survive 
for wfiOO Myrs or longer. Again, clouds of higher column 
densities can have much longer time-scales and can prac¬ 
tically survive photionzation from the cosmic UV back¬ 
ground. This is supported by the fact that some of the 
HCGs show large HI structures within w hich TDG have 
been detected even a distances of 100 kpc (|de Mello et al.l 
l2008f) . Without detailed hydrodynamical simulations it 
is not clear the exact survival times for diffuse gas in each 
of these groups. The conditions such as the temperature 
of the IGrM, the dynamical mass and effective galaxy 
densities within the virial radii, total gas content of the 
IGrM etc. are important parameter in determining the 
lifetime of the diffuse gas structures. However, we be¬ 
lieve that in most groups the basic arguments presented 
in section 3.5.1 will hold and the diffuse gas will survive 
for atleast 500 Myrs. 
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4. CONCLUSIONS 

We presented 21cm GBT H I observations of HCGs 31, 
44, 68, and 92. The observations were carried out in 
pointing map modes resulting in a 5x5 map covering an 
area of 25'x25'. The pointings were separated by 4'. The 
central pointings we re positioned at the center of the 
groups as defined bv iHickson et al.l (Il992f) . 

1. We detected H I in all the four groups. The 
strengths of H 1 emission and the spectral shapes 
observed at the central pointing are consistent with 
our previous single dish measurements. The total 
H 1 masses were found to be 1.86 x 10^°, 1.82 x 
10^, 8.32 X 10®, and 2.45 x lO^® Mq for HCGs 31, 
44, 68, and 92 respectively. 

2. The tidal debris-dominated groups showed en¬ 
hanced blending between various spectral features 
due to the low velocity dispersion among the fea¬ 
tures. Most of the H 1 in these groups exists in 
single peaked Gaussian-like profile center at the 
redshifts of member galaxies. HCG 92 shows high 
velocity wings that connect between two distinct 
spectral features. For ISM-dominated groups, most 
of the H I lies in disks as confirmed by doubled- 
horned profiles (44D and 68C). 

3. We found the H 1 observed with the GBT to have 
the same spatial distribution as the high-surface 
brightness gas seen with the VLA for HCGs 31 and 
68. More than 95% of the total H 1 in these two 
groups is centrally concentrated as indicated by the 
minimal difference in flux between the total H I and 
the pointing with maximum flux. The H 1 distri¬ 
butions in HCGs 92 and 44 were spread out and 
covered most of our maps. The total H 1 in these 
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two groups was 23% and 85% higher than that ob¬ 
served in the pointing with maximum flux. 

4. We modeled expected the H I masses for each 
pointing based on the VLA H I map. While 
HCGs 31 and 68 show very little deviation from 
the modeled values, HCGs 44 and 92 show signifi¬ 
cant offsets. The faint gas fraction for HCG 44 is 
large in the Northwest-Southeast direction below 
the group center. Similarly in HCG 92, we found 
diffuse gas seen in the H I “wing” between veloci¬ 
ties 6300-6530 km s“^ is located in the northwest 
of the group. The detection of the wing feature in 
only a few pointings indicates that it is spatially 
localized and near Arc-N. 

5. The spatial and dynamical similarities between the 
total (faint-l-HSB) and the HSB gas indicate that 
the faint gas detected by the GBT is of tidal ori¬ 
gin, like the HSB component. We investigated the 
effects of the cosmic UV background and the escap¬ 
ing ionizing photons from the star forming regions 
on the survival of the neutral gas. We found that 
the gas will survive and stay primarily neutral for 
at least 500 Myrs but possibly much longer. 
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TABLE 1 

Properties oe our Sample. 


Group 

Redshift 

Zgp 

Velocity 

Up ^ 

(km s 

Distance 

(Mpc) 

Group diameter 

^9P ^9P 

(arcmin) (kpc) 

Map center 

R.A. Dec 

Scale 

(kpc/') 

Previous HI imaging studies 

HCG 31 

0.01347 

4039 

58.3 

4.2“ 

69 

05:01:38.43 -04:15:25.0 

16.5 

Verdes-Montenegro et al. f2005) 

HCG 44 

0.00460 

1379 

19.8 

16.4 

93 

10:18:00.45 -1-21:48:42.9 

5.7 

Williams et al. (19911: Serra et al. (20131 

HCG 68 

0.00800 

2310 

34.5 

9.2 

91 

13:53:40.99 -1-40:19:09.5 

9.8 

Tang et al. (2008) 

HCG 92*^ 

0.02150 

6446 

93.6 

3.2 

84 

22:35:57.52 -1-33:57:37.4 

26.1 

Williams et al. l'2002'l 


^ The group diameter for HCG 31 was revised from the original value of Ogp — 0.9^ lIHicksonll 19821 1 to 9gp — 4.2' by IVerdes-Montenegro et al.l 11200511 based 
on VLA HI observations. 

HCG 92 is known as Stephan’s Quintet. Our observations do not include the foreground galaxy HCG 92A (aslo known as NGC 7320) at a redshift of 
0.002622 (u ^ 786 km s“^ ). 


TABLE 2 

Details of the Observations. 


Group 

integratzon 

(s) 

AV^i 
(km s 1 ) 

HI mass limit'^ 
(Log[M 0 ]) 

HGG 31 

480 

3900-4250 

8.69-8.86 

HGG 44 

260 

1000-1750 

7.88-7.91 

HGG 68 

420 

2100-2500 

8.30-8.46 

HGG 92 

300 

5950-6950 

9.26-9.33 


^ On-source integration per pointing. 

^ The range of velocities over which the HI emission 
was detected. There may be multiple components found 
within this velocity window. 

^ The values show the range of Msg. for all the pointings 
around each of the HCG. The HI mass was estimated as¬ 
suming a line width of AVhi =400 km s“^ in line-free 
spectral windows. 


TABLE 3 

Angular distances (in arcminutes) 

OF POINTINGS FROM THE CENTER OF 
THE MAP. 


1 1 2 3 4 5 


A 

11.3 

8.9 

8.0 

8.9 

11.3 

B 

8.9 

5.7 

4.0 

5.7 

8.9 

C 

8.0 

4.0 

0.0 

4.0 

8.0 

D 

8.9 

5.7 

4.0 

5.7 

8.9 

E 

11.3 

8.9 

8.0 

8.9 

11.3 
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TABLE 4 

HI MASSES FOR EACH POINTING. 


Position 

HCG 31 

(Log[MQ]) 

HCG 44 
(Log[M0]) 

HCG 68 
(Log[MQ]) 

HCG 92 

(Log[MQ]) 

Al 

<8.56 

_ 

<8.30 

<9.27 

A2 

9.06 

- 

8.73 

9.65 

A3 

9.29 

8.58 

9.37 

9.66 

A4 

9.09 

8.64 

9.55 

<9.29 

A5 

<8.63 

8.26 

9.25 

<9.28 

B1 

8.88 

8.25 

<8.35 

9.84 

B2 

9.76 

8.52 

9.19 

10.17 

B3 

10.01 

8.87 

9.75 

10.18 

B4 

9.80 

8.99 

9.91 

9.60 

B5 

9.02 

8.72 

9.64 

<9.27 

Cl 

9.37 

8.35 

<8.35 

10.13 

C2 

10.01 

8.79 

9.09 

10.30 

C3 

10.26 

8.97 

9.62 

10.20 

C4 

10.04 

8.98 

9.80 

9.59 

C5 

9.24 

8.72 

9.58 

<9.26 

D1 

9.25 

<7.91 

<8.39 

9.64 

D2 

9.90 

8.46 

<8.40 

9.68 

D3 

10.11 

8.70 

9.08 

9.77 

D4 

9.83 

8.87 

9.40 

<9.27 

D5 

8.90 

8.50 

9.13 

<9.27 

El 

9.10 

<7.89 

<8.40 

<9.28 

E2 

9.42 

<7.88 

<8.41 

<9.28 

E3 

9.47 

8.25 

<8.45 

<9.29 

E4 

9.04 

8.63 

<8.45 

<9.32 

E5 

<8.84 

8.59 

8.44 

<9.32 

Total“ 

10.27 

9.26 

9.92 

10.39 


Estimated from composite HI spectra shown in Figure 3. 










